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The Schrodinger equation with a two-dimensional 5-function potential is a simple 
example of an asymptotically free theory that undergoes dimensional transmuta- 
tion. Renormalization requires the introduction of a mass scale, which can be 
lowered perturbatively until an infrared cutoff produced by non-perturbative ef- 
fects such as bound state formation is encountered. We outline the effective field 
theory and similarity renormalization group techniques for producing renormalized 
cutoff hamiltonians, and illustrate the control of logarithmic and inverse-power-law 
errors both techniques provide. 



1 Motivation and Outline 

This article follows a more complete discussion in a longer article being pre- 
pared for publication!!! Elsewhere in these Proceedings the possibility of de= 
riving a valid constituent approximation in light-front QCD is discussedfl 
A simple procedure that produces a renormalized light-front QCD hamilto- 
nian starts with the computation irf a cutoff hamiltonian using a similarity 
renormalization group (SRG) autffi and coupling coherence. Lffim In this arti- 
cle we illustrate this step using the two-dimensional <5-function in a one-body 
Schrodinger epuationJiB We also discuss a related approach, effective field 
theory (EFT) OB 

Locality is the source of perturbative divergences in field theory, and it 
produces similar effects in simple one-body quantum mechanics problems. Lo- 
cality leads us to consider potentials consisting of S- functions and their deriva- 
tives. These potentials produce divergences that require the introduction of 
a mass scale. Changing this mass scale should be equivalent to changing a 
dimcnsionless coupling constant (i.e., dimensional transmutation) for a scale- 
invariant local hamiltonian. The two-dimensional <5-function problem is solved 
in the second section, providing the data for our error analysis of EFT and 
SRG calculations. 

Problems requiring renormalization result from an arbitrarily large num- 
ber of degrees of freedom being coupled. This is clear in momentum space, 
where the off-diagonal matrix elements of a 5-function are constants. Diver- 
gences result from an infinite number of scales being directly coupled, and 
renormalization requires a cutoff A to regulate the divergences. Physical re- 
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suits such as binding energies should be independent of A, so there should 
exist a line of cutoff hamiltonians H\ that all produce the same physical re- 
sults. A renormalization group (RG) is built from transformations, T(A, A'), 
that change the cutoff. T connects hamiltonians at different scales and in the 
SRG T is typically a unitary transformation. 

In EFT T exists in principle, but H\ is determined at each scale by con- 
straining the couplings in front of a finite number of local operators to fit 
data. A simple scale transformation leads to the classification of these oper- 
ators as relevant (powers of the cutoff appear), marginal (the cutoff appears 
only in logarithms that are absorbed in the running marginal coupling), and 
irrelevant (inverse powers of the cutoff appear) . The same operators arise in 
the SRG, where their strengths are approximated by expansions in powers 
of a running marginal coupling. For the simple two-dimensional (5-function 
there are no relevant operators and only one marginal operator. There are 
always an infinite number of irrelevant operators, and these become increas- 
ingly important as the cutoff approaches a non-perturbative energy scale in 
the problem. 

The easiest way to clarify these points is by example. We show a complete 
solution of the two-dimensional (5-function problem, introducing the K-matrix 
to compute scattering phase shifts. There is one bound state, and its binding 
energy sets the scale for non-perturbative physics. The EFT hamiltonian is 
produced by truncating the series of irrelevant operators, and then fitting 
a finite number of remaining couplings to an equal number of low-energy 
scattering constraints. This leads to power-law suppression of errors in the 
predicted binding energy, with the expansion clearly breaking down as the 
cutoff approaches the binding energy. The SRG hamiltonian is produced 
by solving a differential equation for H\ subject to the constraint that all 
irrelevant couplings are analytic functions of the single marginal coupling. 

2 The two-dimensional (5-function potential 

Consider the Schrodiager equation in two dimensions with an attractive Dirac 
(5-function potentiality 

- V^(r) - a <5^(r) *(r) = E *(r) . (1) 

The coupling ao is dimensionless, so the hamiltonian is scale invariant (i.e., 
there is no intrinsic length or energy scale). It is relatively straightforward to 
solve this equation in position space by introducing a convenient distribution 
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function to regulate the (^-function, but we will work in momentum space to 
stay close to the EFT and SRG calculations. 

The Schrodinger equation in momentum space is, 

p 2 Hp) - ^3 / A *(q) = E *(p) • (2) 

As a consequence of scale invariance, if there is a bound state solution to Eq. (||) 
then it will admit solutions for any E < 0. This corresponds to a continuum 
of bound states with energies extending down to -co, so the system is not 
bounded from below. By rearranging the terms in the Schrodinger equation 
we obtain 

. a *(0) 

where \l/(0) is the position space wave-function at the origin and Eq = —E > 
is the binding energy. 

To obtain the eigenvalue condition for the binding energy, we can integrate 
both sides of Eq. (§): 



1 = 7^ / *P73-TTrv. (4) 



poo. 1 

a f 1 

^Jo dpp WTE ) 

The integral on the r.h.s. diverges logarithmically, so the problem is ill- 
defined. 

Renormalization requires that we regulate the divergence and absorb cut- 
off dependence in the coupling. First, we regulate the integral with a momen- 
tum cutoff, obtaining 

1 = £* ^ fry 1 -g^fl.^ , (5) 

2tt Jo PP (P 2 + E ) 4vr V Eq) ' W 

so that 

A 2 

E = — ■ (6) 

e a ° — 1 

Clearly, if the coupling cvq is fixed then Eq — > oo as A — ► oo. In order 
to eliminate the divergence and produce a finite, well-defined bound state we 
can renormalize the theory by demanding that the coupling runs with the 
cutoff A in such a way that the binding energy remains fixed as the cutoff is 
removed: 

47T 

"o -> "a = — 7 — v ■ (7) 



hi 



japanrg: submitted to World Scientific on February 9, 2008 3 



The dimensionless renormalized running coupling a\ that characterizes 
the strength of the interaction is therefore replaced by a new (dimensionful) 
parameter Eq, the binding energy of the system. This is a simple example of 
dimensional transmutation: even though the original "bare" hamiltonian is 
scale invariant, renormalization leads to a scale that characterizes the physical 
observables. Note that Eg can be chosen arbitrarily, fixing the energy scale of 
the underlying (renormalized) theory. It is also interesting to note that the 
renormalized running coupling a\ vanishes as A —> oo and so the theory is 
asymptotically free. 

This renormalized hamiltonian can be used to compute other observables. 
The usual prescription for the calculations is to obtain the solutions with the 
cutoff in place and then let it approach oo, using the running coupling in 
Eq. 0. If an exact calculation can be implemented, the final results should 
be independent of the regularization and renormalization schemes. As an 
example, we calculate the scattering wave function, 



$k(p) =^ (2) (P- 



*(0) 



2tt (p 2 — k 2 + i e) 



(8) 



where k — V E. Integrating both sides over p with a cutoff A in place, we 
obtain 



*(0) 



1 

2^ 



A 2 



-k 2 



(9) 



thus, 



«A *(0) = ^ 



1 



J_ l n ( l + - — l r; j | 



A 2 



-k 2 + ie 



(10) 



In the limit A — > oo we obtain: 

a A *(0) 



The resulting scattering wavefunction is then given by 



(11) 



$ fc (p) = ^(p-k) + -L 



2tt (p 2 -k 2 ) [ \Eo 



(12) 



Only S-wave scattering occurs, because the centrifugal barrier completely 
screens the ^-function potential for higher angular momentum. 
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For our numerical calculations we used the K-matrix, which has mixed 
boundary conditions that lead to a principal value. The K-matrix Lippmann- 
Schwinger equation with the renormalized potential is given by: 

K (p, p'; k) = V(p, p') + V J d 2 q tf(q, p'; k) . (13) 

Since only S-wave scattering takes place we can integrate over the angular 
variable, obtaining 

K^(p,p';k) = V^( P y)+r£ dq q K^\q,p';k) , 

(14) 

where 

V^°\p,p') = -^. (15) 
The Lippmann-Schwingcr equation for the "on-shell" K-matrix is given 

by: 

Solving this equation and taking the limit A — > oo, we obtain the exact on- 
shell K-matrix: 

K ^ = -—hx- 



Using 



k cot 6 (k) =-i*_i_, (18) 
7T K {k) 



we can obtain the exact phase-shifts: 



cot S = - In ( — ) . (19) 



3 Effective field theory 

3.1 EFT formalism 

EFT provides a systematic procedure for improving the approximate descrip- 
tion of composite effective degrees of freedom and their dynamics by increasing 
the number of local effective interactions. The effects of high-energy degrees of 
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freedom are incorporated in effective interactions whose strength is adjusted to 
fit appropriate low-energy results. The keys to maintaining predictive power 
in EFT are approximate locality and an expansion in powers of a small ratio 
of low-to-high energy scales or equivalent^ short-to-long distance scales. 

We follow Lepage's treatment of EFT.O For simplicity we use a separable 
hamiltonian with a cutoff A that smoothly regulates interactions involving 
high energy states: 



Vx( P ,p') 



C (A)+C 2 (A) 



(p 2 + P ' 2 ) 
2A 2 



(p 4 +p' 4 ) 
4A 4 



2 /2 

p p 
2A 4 



x e 



(20) 



The parameters Ci in the expansion can be determined by fitting data for 
low-energy processes. Here we use as data the values for the inverse on-shell 
K-matrix, calculated witb_the exact theory. We follow the method described 
by Steele and Furnstahl 113 with the parameters fixed so that EFT gives the 
same momentum expansion for the inverse on-shell K-matrix as the exact 
theory to a given order. 

The Lippmann-Schwinger equation for the K-matrix with is 



K x (p,p';k)^V x (p,p , )+V / d 2 q 



Vx(p,q) 
k 2 -a 2 



Kx(q,p';k) 



(21) 



where V denotes the principal value. The equation refers to the S-wave poten- 
tial (angular variable integrated out) and d 2 q = q dq. To obtain IJm-K-matrix 
we solve the Lippmann-Schwinger equation nonperturbativelyO'tjLS The ef- 
fective potential given in Eq. (f20|) (truncated at order p A ) can be written in 
the form 



V\{jp,p') = e 



X> 2i 

i,j=0 



A 



ij P 



>2j 



(22) 



where Ay are the matrix elements of 



/ C (A) C 2 (A)/2A 2 C 4 (A)/4A 4 \ 



C 2 (A)/2A 2 C^(A)/2A 4 
\C 4 (A)/4A 4 







(23) 



/ 



The solution of the Lippmann-Schwinger equation can then be written in the 
form 
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K\(p',p; k) = e e fer ^ p 21 r# (k) p 

i,3=0 

The unknown matrix r(fc) satisfies the equation 

r(k) = A + A J(fe) r(k) , 

where 

( h{k) h(k) I 2 {k)\ 



h(k) h{k) h(k) 

\l 2 {k) Jg(fc) / 4 (fc)/ 



(24) 



(25) 



(26) 



with 



In(k) 



dq q 



2n+D-l 



k 2 



(27) 



Solving Eq. ( p5| ) for r analytically and substituting the solution in Eq. ( p4| ) 
we obtain the effective K-matrix. 

It is straightforward to extend the calculation to higher orders. We fit 
the difference between the effective and exact inverse on-shell K-matrix (p 2 = 
p' 2 = k 2 ) to an interpolating polynomial in k 2 /X 2 to the highest possible 
order, 



A 



= A + A s 



A 2 



Ai 



A 4 



(28) 



The coefficients Ai are minimized with respect to the variation in the param- 
eters Ci of the effective potential. The number of coefficients that can be 
minimized is given by the number of parameters appearing in the effective 
potential. 

By adjusting only the parameter Co (A) we should eliminate the leading 
error. As each term is added to the effective potential, followed by the ad- 
justment of the respective parameter, we expect the errors in the on-shell 
K-matrix to be systematically reduced by powers of k 2 /A 2 . Thus, in a log-log 
plot for A[l/K] we expect to obtain straight lines with slope given by the 
dominant power of k/X in the error. 

The expansion given by Eq.( |28| ) becomes invalid at best when the mo- 
menta involved are of the same order as the cutoff, where one expects the short 
distance effects to be directly resolved. This point corresponds to the radius 
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of convergence of the effective theory. Using the effective potential (truncated 
at a given order) with the parameters fixed by fitting the inverse on-shell 
K-matrix we diagonalize the effective hamiltonian numerically, obtaining a 
bound-state energy that can be compared to the exact result. 



3.2 EFT results for the two-dimension 5 -function 

The EFT power counting scheme is simple. For the leading-order prediction, 
we use the effective potential with one parameter. Evaluating the difference 
A [l/i^ )] = l/K^' — 1/Kq and expanding it in powers of k/X we obtain: 



1 




' 1 


7 


_X(0)_ 




.Co 


2 



By choosing 



hi 



A^ 
E 



1 

Co 



111 

A 2 



O 



(29) 



Co (A) = 



In 



(*) 



(30) 



7 



the leading logarithmic error is eliminated and the remaining errors are dom- 
inated by the term 



O 



7 — In 



A^ 

E 



1 



A 2 



(31) 



To evaluate the errors in the binding energy we use the effective potential 
with one, two and four parameters adjusted to fit the inverse K-matrix. We 
solve the Schrodinger equation numerically for different values of the cutoff 
A. In Fig. 1 we show in a log-log plot the absolute values for the relative 
errors in the bound-state energy as a function of Eo/X 2 . As expected, we 
see straight lines with slope given by the dominant power of Eg /A 2 in the 
error. Note that the dominant logarithmic errors are completely absorbed 
by adjusting the parameter Co and when more parameters are adjusted we 
obtain power-law improvement. 



4 The similarity renormalization group 

4-1 SRG formalism 

The similarity renormalization group i@ is made of similarity transformations 
that run a smooth cutoff on the magnitude of off-diagonal matrix elements. 
As the cutoff is lowered, the hamiltonian is forced toward the diagonal. It 
is straightforward to show that an infinitesimal transformation produces a 
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differential equation governing the evolution of the hamiltonian. Glazek and 
Wilson considered very general transformations in their original work,B but 
we use the simple transformation developed by Wegnerfl The infinitesimal 
transformation is 

^ = [H s ,[H s ,h]. (32) 

Here H s = h + V s , where h is the free hamiltonian, evolves with a flow 
parameter s that ranges from to oo. The flow-parameter has dimensions 
1/ (energy) 2 and can be expressed in terms of the similarity cutoff A by the 
relation s — 1/A 2 . Wegner advocates using the full diagonal part of the 
hamiltonian instead of using h in the transformation. 

The reduced interaction, V S ij (the interaction with a gaussian similarity 
cutoff factored out) is defined by 

V sij = e" sA -- V sij , (33) 

where Ay = p| — p 2 for our simple transformation. The gaussian factor forces 
the interaction toward the diagonal to leading order, and the equation for the 
evolution of V starts at second order. 
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The flow equation for the reduced interaction is 



si] 



+ A jk ) V slk V skj e - 2sA ^ , (34) 



ds 

k 

where we use - Af k - A 2 k = -2A ik A jk . 

To solve this equation we impose a boundary condition, H s \ g ^ g = H So , 
where so can be expressed in terms of the large similarity cutoff as so = 1/A 2 . 
Next we make a perturbative expansion, 

_ _(1) _(2) 

V S = V S +V S + • • • , (35) 

where the superscript implies the order in the running canonical coupling. All 
counterterms are determined using coupling coherence, which assumes V is a 
perturbative function of a. 



At first order we have 



which implies 



dV... 

-^=0, (36) 



_(!) 



V sij = Vs oij , (37) 

where s is the final scale. Because of the dimensions of the flow parameter 
we have s > s , corresponding to a smaller cutoff. The "no cutoff limit" 
corresponds to s — ► 0. 

At second order we have 

_(2) 

=E + V ^ V °oV e- 2sA > kA > k • (38) 

k 

Integrating, we obtain 

— < 2 > _ / 1 1 

The rcnormalizcd hamiltonian for the non-relativistic delta-function po- 
tential is given by 



H x (p, p') = p 2 S^(p p') + e- {J ^ [vfW) + V^(p, p') + ... 



(40) 
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where 



vf } (p.p') 



a\,i 



(41) 



vfW) 

^ (n) (p,p') 



(2tt)* ' 

«L^ (2) (p,p'), 
^f. w (p.p')- 



(42) 
(43) 



Here A is a momentum cutoff related to the flow parameter by s = 1/ A and 
i denotes the order of the calculation for the running coupling. 

The renormalized hamiltonian can be used to compute eigenvalues and 
eigenstates. Since the hamiltonian is derived perturbatively we expect cutoff 
dependent errors in the observables. Formally, we can regroup the terms in 
the renormalized hamiltonian and write it as a momentum expansion similar 
to the one for the EFT in Eq. (pp|). The difference is that the expansion 
parameters are perturbative functions of the running coupling a\. 

There are two interdependent sources of errors in the perturbative sim- 
ilarity renormalization group. First, a\ is approximated and this leads to 
inverse-logarithmic errors. Second, the coefficients of all irrelevant opera- 
tors are approximated perturbatively, and this leads to products of inverse- 
logarithmic and inverse-power-law errors. 

In our examples we focus on the bound state errors. We fix the coupling 
at one scale and use the flow-equation to obtain the coupling as a function 
of the cutoff A to a given order. We then perform a sequence of bound- 
state calculations with better approximations for the hamiltonian. Once the 
sources of errors are identified, it becomes relatively simple to analyze order- 
by-order how such errors scale with A. In principle, to completely eliminate the 
errors proportional to some power m in the momentum expansion we should 
use the similarity hamiltonian with the exact running coupling (renormalized 
to all orders) and include the contributions up to 0(p m / X m ) coming from 
all effective interactions (all orders in ax)- We should emphasize that in 
a realistic calculation we would fit the coupling a\ to an observable. This 
nonperturbative renormalization eliminates the dominant source of errors we 
display in SRG calculations in this article. We choose to renormalize the 
coupling perturbatively here because the only observable we compute is the 
single bound state energy of a delta-function potential, and fitting this energy 
would prevent us from displaying errors. 

4-2 SRG results for two-dimensional d-function 

Glazek and Wilson have studied a discretized analog of the two-dimensional 
6 function, which allowed them to perform both perturbative and nonpertur- 
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bative transformation^ Their perturbative results are similar to ours. In the 
two-dimensional case the canonical hamiltonian in momentum space with a 
delta-function potential can be written as 

H(p,p')=h(p,p')+V(p,p'), (44) 

where h(p, p') — p 2 S^(p — p') corresponds to the free hamiltonian and 
y(P)P') = — a o/(27r) 2 corresponds to the Fourier transform of the delta- 
function potential. 

Integrating out the angular variable, the flow equation obtained with 
Wegner's transformation in terms of matrix elements in the basis of free states 
is given by 

dVs ^ P>) = -(P 2 ~P' 2 ) 2 V s (p,p')- r dk k (2fc 2 -p 2 "P' 2 ) V s (p,k) V s (k,p') . 
ds J 

(45) 

In principle, we can set the boundary condition at s = (no cutoff), i.e, 

H s=Q (p,p') = H(p,p')=p 2 6^(p- P ') - g . (46) 

However, the hamiltonian with no cutoff produces logarithmic divergences. 
The boundary condition must be imposed at some other point, leading to 
dimensional transmutation. The reduced interaction V s (p,p') is defined such 
that 

V s (p,p') = e- S ^-^ 2 V s (p,p>) . (47) 

Assuming that h is cutoff independent we obtain the flow equation for the 
reduced interaction, 

dZ = _^ 2s pV 2 r dk k {2k 2 - p 2 - p' 2 ) e -M^-^(^+ P ' 2 )] 
as J 

xV s (p,k)V s (k,p') . (48) 
This equation is solved using a perturbative expansion, starting with 

K (1 W) = "g- (49) 

We assume a coupling-coherent solution in the form of an expansion in powers 
of a s /2-K, satisfying the constraint that the operators Fg n \p,p') vanish when 
p = p' = 0, 

OO 

K(p,p') = -g + E(£)" ^i"W)- (50) 



japanrg: submitted to World Scientific on February 9, 2008 12 



Note that the expansion parameter is a s /27r. 

Using the solution Eq.(j50|) in Eq.(Es) we obtain 



dV s 1 da s 

~ds~ = ~ (2tt) 2 ~~ds~ 



1 



oo 

n=2 v ' 
I2\ 



-1 ^1 F (n) (j) v n , a n dF^\p,p') 
ds s [p > p) + a ° ds 



dk k [2k z -p z -p' z ) e 



-2 S bV 2 +fc 4 -fc 2 (p 2 +p' 2 )] 



_i_ ( as 



2n 



F^(p,k) 



(51) 



This equation is solved iteratively order-by-order in a s /2n. Again, if a s /2ir 
is small the operator V^ 1 (p,p') can be identified as the dominant term in the 
expansion of V s (p,p') in powers of p and p'. This operator corresponds to a 
marginal operator (since the coupling is dimensionless and there is no implicit 
mass scale). The higher-order terms correspond to irrelevant operators. 

For simplicity we use the interaction with a large cutoff on all momenta, 
A, to define the exact theory. We define 

Air 

a So=0 = QA = — t — y > (52) 



ml 



A 2 . 

So 

and set all irrelevant operators to zero at s = 0. Note that the coupling a\ Q 
is fixed at Ao = oo by fitting the exact binding energy. With this definition 
the similarity hamiltonian with no cutoff becomes well-defined and we can 
set all of the similarity transformation boundary conditions at so = 0. All 
integrals below are cut off (p <C A) 
At second-order we have 



1 ckt s ,A 
2?r ds 



1 



(2ny 



ds 



(53) 



where 



i 



(27T) 2 
1 



dk k (2k 2 

2s p 2 p' 2 



p^e- 



ls 



(54) 



The resulting second-order running coupling and irrelevant operator are given 
respectively by 

OA 



a s ,A.2 



1 - f£ [7 + In (2sA 4 ) - Ei (-2sA 4 )] 



(55) 
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Figure 2. The SRG bound-state energy using various approximations for the similarity 
hamiltonian. The exact theory is fixed by choosing A = 50 and Eq = f . 



and 

*1 2 aW) = \ b + ln (2* P 2 P' 2 ) - Ei(-2 S p 2 p' 2 )] 
+ + ln(2sA 4 ) - Ei(-2sA 4 )] 

- \ [ 7 + In (a [(p 2 A 2 ) 2 + (p> 2 - A 2 ) 2 - (p 2 - p' 2 ) 2 } ) 

-Ei (s [(p 2 - A 2 ) 2 + (p' 2 A 2 ) 2 (p 2 -p 12 ) 2 ])] . (56) 

In Fig. 2 we show the binding-energy as a function of the cutoff and in 
Fig. 3 we show the errors in the binding-energy obtained using the following 
approximations for the potential with A = 50: 

(a) marginal operator with coupling (a\ ), 

V x ,^v') = -^e-^^ ; (57) 

(b) marginal operator with running coupling renormalized to second-order 

t r i m C*A.A,2 _(p 2 -p' 2 ) 2 

Vx, A (p,p ) = J- 2 - e 13 ; (58) 
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Figure 3. The SRG errors in the binding energy using various approximations for the sim- 
ilarity hamiltonian. The exact theory is fixed by choosing A = 50 and Eq = f . 



(c) marginal operator plus second-order irrelevant operator with running cou- 

(2) 

pling renormalized to second-order (aA,A,2i F\ )j 



Vx,a(p,p') 



Ot\\.2 



,A,2 . /aiA,2\ z?(2)/ /\ 



i 2 /2 

(p -v 



(59) 



(d) marginal operator with running coupling renormalized to third-order 

( Q A,A,3), 



— — e 

2tt 



(60) 



(e) marginal operator plus second-order irrelevant operator with running cou- 

(21 

pling renormalized to third-order ax,A,3) F\ )> 



^,a(p,p') = 



QA,A,3 

2tt 



«A,A,3 

2tt 



/ 2 /2,2 
(P — P ) 



(61) 



In Fig. 2 we clearly see all of the binding energies approaching the exact 
result, which we choose to be 1. If the coupling does not run (a) the energy falls 
to zero as the cutoff is lowered. When the coupling is allowed to run to second 
order (b) the binding energy deviates from the exact result more slowly as the 
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cutoff is lowered, and subsequent improved approximations further remove 
dependence on A. 

Fig. 3 is a log-log plot of the relative error in the binding energy as 
a function of Eo/X 2 . When A 3> A the errors scale like inverse powers of 
E n /X 2 . They are forced to approach the correct result as A — > oo, and the 
exact running coupling deviates from the approximate running coupling as 
a power of A/A in this region, even when the approximate coupling is not 
allowed to run. The coefficients of these powers are large, and the analysis 
changes as A — > A. The errors become logarithmic for (a), but they scale 
as an inversedogarithm squared in (b) where the coupling runs to second- 
order. This removes the dominant inversedogarithmic errors, and we can 
clearly see further improvement when the coupling runs to third-order in (d). 
The addition of irrelevant operators in (c) and (e) are important, but the 
third-order error in the running coupling becomes as large as the error from 
the leading irrelevant operators, which shows up when (d) crosses (c). When 
the second-order irrelevant operators are added with the coupling running to 
third order, we obtain approximate power-law improvement over (b). 

5 Conclusion 

Renormalization not only removes unphysical divergences, it also allows us 
to work with finite cutoff's while controlling the errors they introduce. Finite 
cutoffs are often required for nonperturbative calculations, and the precise 
control of the errors they introduce is essential. 

In light-front QCD we can use the SRG to compute the hamiltonian order- 
by-order in the canonical QCD coupling, a\. We are free to choose A so 
that multi-parton high-energy states effectively decouple from few-parton low- 
energy states, as long as this happens before a\ becomes too large. We adjust 
the coupling non-perturbatively by fitting low-energy data. This should re- 
move all inverse-logarithmic errors, leaving only powers of inverse-logarithms 
and inverse-powers of the cutoff. We can systematically lower the errors by 
adding irrelevant operators, all of which are fixed by coupling coherence. 

Acknowledgments 

One of us (R.J.P.) is very grateful to the hosts of the 1998 YITP- Workshop 
on QCD and Hadron Physics. The opportunity to see Kyoto and to interact 
with so many excellent young Japanese physicists was fantastic. We would also 
like to acknowledge many useful discussions with Brent Allen, Dick Furnstahl, 
Stan Glazek, Roger Kylin, Rick Mohr, Jim Steele, and Ken Wilson. This work 



japanrg: submitted to World Scientific on February 9, 2008 16 



was supported by National Science Foundation grant PHY-9800964, and S.S. 
is a CNPq-Brazil fellow (proc. 204790/88-3). 

References 



10. 

11. 
12. 



13. 



Sergio Szpigel and Robert J. Perry, in preparation. 
Robert J. Perry, these Proceedings (1999). 
S. D. Glazek and K.G. Wilson, Phys. Rev. D 48, 5863 (1993); 
D. Glazek and K.G. Wilson, Phys. Rev. D 49, 4214 (1994). 
Wegner, Ann. Physik 3, 77 (1994). 

D. Glazek a nd K.G. Wilson Phys. Rev. D 57, 3558 (1998). 
S. Walhout, |hcp-th/98060"97l (1998). 

Oehme and W. Zimmermann, Commun. Math. Phys. 97, 569 (1985); 



Oehme, K. Sibold and W. Zimmermann, Phys. Lett. B 147, 115 
(1984). 

R. J. Perry and K. G. Wilson, Nucl. Phys. B 403, 587 (1993); 
R. J. Perry Ann. Phys. 232, 116 (1994). 

B. H. Allen and R. J. Perry Phys. Rev. D 58, 125017 (1998), |hcp- 



th/9804136 



B. H. Allen, 

in M. A. B. Beg Memorial Volume, A. Ali and P. Hoodbhoy, 



and R. J. Perry, |hep-th/9812080 
al \ 
1991) 



R. D. Kylin 
R. Jackiw 

eds. (World S cientific, Singapo re 

G. P. Lepage, |aucl-th/9706029 . 

S. Weinberg, Physica A 96, 327 (1979); 

H. Gcorgi, Ann. Rev. Nucl. Part . Sci. 43, 209 (1993); 
V. Mano har, jhc-p-ph/9606222] ; 

. B. Wise, |hep-ph/9707522|; 
B. Kaplan, |nucl-th/950603 



A. 
M. 
D. 
U. 
J. 



van Kolck, [nucl-th/9808007 
V. Steele and R. J. Furnstahl, 



Nucl. Phys. A637,46 (1998), nucl- 



th/9802069t 

V. Steele, [nucl-th/9807022] . 
R. Phillips, S. R. Beane and 



J 

14. D 

(1998), |hep-th/97060~70|. 

15. M. C. Birsc , |nucl-th/9804028| ; |nucl-th/9708435 



T. D. Cohen, Ann. Phys. 263, 255 



ph/980839S 



hcp-ph/9807302| , |hep- 



16. J. Gegelia, |nucl-th/9805006 . 



japanrg: submitted to World Scientific on February 9, 2008 17 



